Jerk, snap, and the cosmological equation of state 



Matt Visser 

School of Mathematical and Computing Sciences, Victoria University of Wellington, 
PO Box 600, Welhngton, New Zealand 

E-mail: matt . visserSmcs . vuw . ac. nz 



p = Po + Ko (p - Po) + 77 — 



[p - p^f + 0[{p - p,f 



Abstract. Taylor expanding the cosmological equation of state around the current 
epoch 

1 d^p 

2 d^ 

is the simplest model one can consider that does not make any a priori restrictions 
on the nature of the cosmological fluid. Most popular cosmological models attempt 
to be "predictive" , in the sense that once some a priori equation of state is chosen 
the Friedmann equations are used to determine the evolution of the FRW scale factor 
a(t). In contrast, a "retrodictive" approach might usefully take observational data 
concerning the scale factor, and use the Friedmann equations to infer an observed 
cosmological equation of state. In particular, the value and derivatives of the scale 
factor determined at the current epoch place constraints on the value and derivatives 
of the cosmological equation of state at the current epoch. Determining the first three 
Taylor coefficients of the equation of state at the current epoch requires a measurement 
of the deceleration, jerk, and snap — the second, third, and fourth derivatives of the 
scale factor with respect to time. Higher-order Taylor coefficients in the equation of 
state are related to higher-order time derivatives of the scale factor. Since the jerk and 
snap are rather difficult to measure, being related to the third and fourth terms in the 
Taylor series expansion of the Hubble law, it becomes clear why direct observational 
constraints on the cosmological equation of state are so relatively weak; and are likely 
to remain weak for the foreseeable future. 
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1. Introduction 



This article develops a "phenomenological" approach to the equation of state [EOS] 
of the cosmological fluid, and investigates what would have to be done in order to 
observationally determine the EOS. Even at the hnearized level, where 

P = Po + 1^0 {p - Po) + 0[{p - pq)\ (1) 

the first nontrivial coefficient in the EOS will be seen to be related to the cosmological 
jerk — the third derivative of the scale factor with respect to time, and thence to 
the third-Older term in the Taylor series expansion of the Hubble law. This is the 
fundamental reason why observational determinations of the EOS are relatively poor, 
and why it is possible to choose so many wildly differing a priori models for the EOS 
that nevertheless give good agreement with the coarse features of the present epoch. 

More generally, if we describe the cosmological equation of state in terms of a Taylor 
series expansion around the current epoch 

El d"p 
-r 



P = PO 



^ n\ dp" 

n=l 

then the n-th order Taylor coefficient 
d"p 



(p-por+o[(p-po) 



N^ 



(2) 



dp" 

will be shown to depend on the (n+2)-th time derivative of the scale factor 
d"+2a(t) 



di"+2 



(3) 



(4) 



term, in the Taylor expansion of the 



+ r'^ {d.Q'^ + sin=^ Q d(\?) 



and thence to the (n+2)-th term, the 0(2;""'"^) 
Hubble law. 

In most attempts at cosmological model building one takes a FRW cosmology 

ds^ = dt^ + a(t)^ . 

^ ' [1 - kr' 

plus the conservation of stress-energy 

pa^ + 3[p + p] (j?a = 0; 
and chooses some a priori equation of state 

P = p{p); or P = p{p); 



(5) 



(6) 



(7) 



to derive p(a), and equivalently p(a). The Einstein equations then reduce to the single 
Friedmann equation, which can be written in the form 



-k, 



(8) 



and used to determine a(t). (See, for example, any standard text such as [1, 2, 3].) 

In contrast, let us assume we have a FRW universe with good observational data on 
a{t) — in Weinberg's terminology we have a good "cosmography" [1]. In this situation 
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we can use the Einstein equations in reverse to calculate the energy density p{t) and 
pressure p{t) via 

^2 kc^- 



SttGjv p{t) = 3c' 



SttGn p{t) = -c 



a kc a 
- + — + 2- 



(9) 



(10) 



Under mild conditions on the existence and nonzero value of appropriate derivatives we 
can appeal to the inverse function theorem to assert the existence of a t{p) or t(p) and 
hence, in principle, deduce an observational equation of state 

p{p) = p{t = t{p)); p{p) = p{t = t{p)). (11) 

In view of the many controversies currently surrounding the cosmological equation of 
state, and the large number of speculative models presently being considered, such an 
observationally driven reconstruction is of interest in its own right. 

Now in observational cosmology we do not have direct access to a{t) over the 
entire history of the universe — we do however have access [however imprecise] to 
the current value of the scale factor and its derivatives, as encoded in the Hubble 
parameter, deceleration parameter, etc. This more limited information can still be used 
to extract useful information about the cosmological equation of state, in particular it 
yields information about the present value of the w-parameter and the slope parameter 
Ko defined as 

dp 



P 

Wo = - 
P 



dp 



(12) 



The value of the w-parameter in particular has recently become the center of 
considerable interest, driven by speculation that < — 1 is compatible with present 
observations. Such a value of wo would correspond to present-day classical and 
cosmologically significant violations of the null energy condition. The associated 
"phantom matter" {almost identical to the notion of "exotic matter" in the sense of 
Morris and Thorne [4]) leads to a cosmological energy density that is future increasing 
rather than future decreasing. (Sec, for example, [5, 6]). If w(t) subsequently remains 
strictly less than —1, this will lead to a "big rip" [7] — the catastrophic infinite expansion 
of the universe in finite elapsed time. 

Unfortunately it is very difficult to measure wo and kq with any accuracy — I 
will make this point explicit by relating the measurement of wq to the deceleration 
parameter, and the measurement of to the "jerk" of the cosmological scale factor — 
the third derivative with respect to time. 

For related comments see references [8, 9, 10, 11, 12, 13]. The "cubic" term 
of Chiba and Nakamura [8] is identical to the jerk, as is the "statefinder" variable 
called r by Sahni et al. [9, 10, 11]. The other "statefinder" variable (called s, not 
to be confused with the snap) is a particular linear combination of the jerk and 
deceleration parameters. Padmanabhan and Choudhury [12] have also emphasised the 
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need for constructing models for the cosmological fluid that are unprejudiced by a priori 
theoretical assumptions. A good recent survey of the status of the cosmological fluid 
is [13]. 



2. Hubble, deceleration, jerk, and snap parameters 



It is standard terminology in mechanics that the first four time derivatives of position 
are referred to as velocity, acceleration, jerk and snap. J In a cosmological setting this 
makes it appropriate to define Hubble, deceleration, jerk, and snap parameters as 



H{t) 



1 da 
a di ' 



1 d^a 
a dt2 



1 d^a 



s{t)^ 



1 d^a 
a dt^ 



1 da 
a di 

1 da' 
a di 

1 da 
a dt 



(13) 
(14) 

(15) 

(16) 



The deceleration, jerk, and snap parameters are dimensionless, and we can write 
a{t) = ao 1 1 + //o (^ - ^o) - ^ go ^o - ^o)' + ^ Jo ^o - ^o)' 



+^ soHt{t-toY + 0{[t-tof)Y 



In particular, at arbitrary time t 



2q) + fccVa^ 



2q) + kcy{H^a^) 



(17) 



(18) 



3(//2 + fccVa2) 3[1 + key (H^a^)] ' 

While observation is currently not good enough to distinguish between the three cases 
k — —1/0/ + 1 with any degree of certainty, there is nevertheless widespread agreement 
that at the present epoch Hoao/c » 1 (equivalent to \flo — 1\ <^ 1) . 
Warning: Prom a theoretical perspective, Hoao/c » 1 is a generic prediction of 
inflationary cosmology — this is not the same as saying that cosmological inflation 
predicts k = 0. What generic cosmological inflation predicts is the weaker statement 
that for all practical purposes the present day universe is indistinguishable from a. k = 
spatially flat universe. If our universe happens to be a topologically trivial k = FRW 
cosmology, then we will never be able to prove it. Simply as a matter of formal logic, all 
we will ever be able to do is to place increasingly stringent lower bounds on Hq qq, but 
this will never rigorously permit us to conclude that A; = 0. The fundamental reason for 

I Jerk [the third time derivative] is also sometimes referred to as jolt. Less common alternative 
terminologies are pulse, impulse, bounce, surge, shock, and super-acceleration. Snap [the fourth time 
derivative] is also sometimes called jounce. The fifth and sixth time derivatives are sometimes somewhat 
facetiously referred to as crackle and pop. 
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this often overlooked but trivial observation is that a topologically trivial k — Q FRW 
universe can be mimicked to arbitrary accuracy by a = ±1 FRW universe provided 
the scale factor is big enough. § In contrast if the true state of affairs is /c = ±1, then 
with good enough data on H^aQ we will in principle be able to determine upper bounds 
which (at some appropriate level of statistical uncertainty) demonstrate that k ^ Q. 
Also note that even in inflationary cosmologies it is not true that H{t)a{t)/c S> 1 at all 
times, and in particular this inequality may be violated (and often is violated) in the 
pre-inflationary epoch. □ 
Now the w-parameter in cosmology is related to the Morris-Thorne exoticity 
parameter [4] which was introduced by them to characterize the presence of "exotic 
matter" , matter violating the null energy condition [NEC] : 

e = ^ = -gn(p) [l + w] = - sign(p) , ^ (19) 

Thus if w < —1 and p > we have ^ < and the NEC is violated. In contrast, if 
w < —1 but p < we have > 0, the NEC is satisfied but the weak energy condition 
[WEC] is violated. That is, "phantom matter" (matter with w < 1) is not quite the 
same as "exotic matter" (for which ^ < 0), but the two are intimately related. 
Accepting the approximation that HoaQ /c ^ 1 we have 

'^o~^T^^o>0; wo^-^^—:^; and ^0-^(1 + ^0); (20) 

so that in this situation the wo-parameter and exoticity parameter ^0 are intimately 
related to the deceleration parameter qg. In particular if u^o < — 1, so that the universe 
is at the current epoch dominated by "phantom matter", we also (because in this 
approximation po is guaranteed to be positive) have ,^0 < so that at the current epoch 
this phantom matter is also "exotic matter". Exotic matter is powerful stuff: Apart 
from possibly destroying the universe in a future "big rip" singularity [7], if the exotic 
matter clumps to any extent there is real risk of even more seriously bizarre behaviour — 
everything from violations of the positive mass condition (that is, objects with negative 
asymptotic mass), through traversable wormholes, to time warps [4, 14, 15, 16, 17]. 

3. Taylor series equation of state 

Linearize the cosmological EOS around the present epoch as 

P = Po + «o (p - Po) + 0[{p- Po)'] ■ (21) 
To calculate we use 

§ If the universe has nontrivial spatial topology there is a possibility of using the compactification scale, 
which might be (but does not have to be) much smaller than the scale factor, to indirectly distinguish 
between k = -1/0/ + 1. 
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where numerator and denominator can be obtained by differentiating tlie Friedmann 
equations for p{t) and pit). It is easy to see tfiat at all times, simply from the definition 
of deceleration and jerk parameters, we have 



SttG 



N 



SttG 



N 



dp 
dt 

dp 
dt 



-6c^H 



kc^ 



2c^H 



kc^ 



leading to 



1 \ l-j^ + k^/{Hlal) 
3 [l + q^ + kcy{Hlal) 

which approximates (using HqUo/c ^ 1) to 



1 

3 



Jo 



(23) 



(24) 



(25) 



(26) 



.1 + go. 

The key observation here is that to obtain the linearized equation of state you need 
significantly more information than the deceleration parameter go! you also need to 
measure the jerk parameter jo- If the only observations you have are measurements of the 
deceleration parameter then you can of course determine Wq = Po/po, but this is not an 
equation of state for the cosmological fluid. Determining wq merely provides information 
about the present-day value of p/ p but makes no prediction as to what this ratio will do in 
the future — not even in the near future. (This point is also forcefully made in [12].) For 
this reason there have been several attempts to observationally determine w{z), the value 
of w as a function of redshift. See for example [12] and [9, 10, 11]. Since 2; is a function 
of lookback time D/c, this is ultimately equivalent to determining w{t) — p{t)/p{t), 
and implicitly equivalent to reconstructing a phenomenological equation of state p{p). 
I prefer to phrase the discussion directly in terms of the EOS as that will make it clear 
what parameters have to be physically measured. In terms of the history of the scale 
factor a{t), it is only when one goes to third order by including the jerk parameter jo 
that one obtains even a linearized equation of state. 

Going one step higher in the expansion, by using the chain rule and the implicit 
function theorem it is easy to see that 

d^ ^ p-Kp 

dfP {pf ^ > 

More generally d"p/dp" contains a term linear in d"p/di". Using the Friedmann 
equations then imphes that d"p/dp" contains a term hnear in d"+^a/di"+^. Specifically 
for the first nonlinear term it is relatively straightforward take explicit time derivatives 
and so to verify that 

d^p il + kc^ l\H?:al]) f . ^ . . 9. , 

= - f, n\ 2 11^2 2^^A So 1 + 90 + Jo 1 + Jo + Mo + % + 9o 1 + 2go 

dp2 Q 6po(l + 9o + kc^/[m%\r ^ 

+ (sq + jo + go + ?oJo) jTT^ \- (28) 
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In the approximation HQao/c ^ 1 this reduces to 

go(l + go) + Jo(l + Jo + 4go + gg) + go(l + 2go) 
dp^o 6po(l + ?o)=^ ■ ^ ' 

As expected, this second derivative depends hnearly on the snap. Higher order 
coefficients can certainly be computed but are increasingly complicated and less 
transparent in their physical interpretation. (Calculations are impractical without the 
use of some symbolic manipulation package such as Maple.) To now make the connection 
between the the Taylor coefficients of the cosmological EOS and the various parameters 
appearing in the Hubble law we will need to likewise perform a similar Taylor expansion 
of the Hubble law. 



4. Hubble law to fourth order in redshift 

Note that this entire section is independent of the use of the Friedmann equations and 
depends only on the use of a FRW geometry. 

The physical distance travelled by a photon that is emitted at time and absorbed 
at the current epoch to is 



D^cJ dt^c{to-t^). 



(30) 



(31) 



In terms of this physical distance the Hubble law is exact 

a{t*) a{to — D/c)' 

but impractical. A more useful result is obtained by performing a fourth-order Taylor 
series expansion, 

«(fo) _ , , HoD 2 + go H^D^ Q{1 + *)) + Jo H^D^ 



a{to - D/c) 



+ 



24 



c 2 6 

,so + 8jo + 36^70 + 6g2 H^D* 



+ 



HnD 



24 

followed by reversion of the resulting series z{D) D{z) to obtain: 



Ho 



1 + 



-2 + 



i + go + 



Jo 
6 



(32) 



(33) 



3 , 5 o 1 . 5 . So 

1 + 2^0(1 + qo) + g% - ^Jo - Y^QoJo - ^ 



+ Oiz") 



This simple calculation is enough to demonstrate that the jerk shows up at third order in 
the Hubble law, and the snap at fourth order. Generally, the 0{z'") term in this version 
of the Hubble law will depend on the n-th time derivative of the scale factor. (Also 
note that one of the virtues of this version of the Hubble law is that it is completely 
independent of k, the sign of space curvature.) 
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Unfortunately physical distance D is typically not the variable in terms of which 
the Hubble law is observationally presented. That role is more typically played by the 
"luminosity distance", di- For instance, Weinberg defines [1] 

L 



(energy flux) 



47r dl' 



(34) 



Let the photon be emitted at r-coordinate r = at time t*, and absorbed at r-coordinate 
r = To at time t^. Then it is a purely geometrical result that 

.2 ^0 _ oo 



dL = a(io) 



(35) 



a{t^) a(to — D/c) 

Thus to calculate diiD) we need ro{D). Recall that for a null geodesic in a FRW 
universe 

"*» c dt 



But 



f{ro) = 



(I/- 



firo). 



sin ^ To 



To k^O; 
sinh~^ To k — —1; 



therefore 



*° cdt 
t,{D) a{t) 



cdt 



(36) 



(37) 



(38) 



To be exphcit 



ro{D) 



sm 



*° cdt 
to-D/c o(^) 

*° cdt 
to-D/c a{t) 



k = +1; 



A; = 0; 



(39) 



sinh 



to 



cdt 



k 



-1. 



I to-D/c 

We can now Taylor series expand for "short" distances. First note that 



ro{D) 



*° c dt 

to-D/c 



k 
3! 



/■*» cdt 

' to-D/c Q'(t) 



+ 



*° c dt 

to-D/c «(^) 



(40) 



and now expand the integral to third order. (We can check, a posteriori^ that this 
retains sufficient accuracy in the dj^ ^ D conversion for determining the Hubble law to 
fourth order.) Then 

l + //o(io-t)+ ^-^Hl {to-tf 



*« cdt 

to-D/c «(^) 



cdt 



to-D/c 

^ 6(1 + go) + Jo ^3 



(to-t)' + C'[(to-t)' 



(41) 
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c 



1 



-{D/c+^Ho{D/cf + 



6(1 + go) + jo 



24 



ao 
+ 

ao [ 2 c 
+0 



2 + go 



(D/c)^ 



(42) 



2 + go 



+ 



6(1 + go) + Jo 

24 



HnD 



(43) 



The conversion from physical distance travelled to r coordinate traversed is given by 



r-oiD) 



D 
ao 



IHqD 1 

1 + — + - 

2 c 6 



2 + go 



1 

+ 24 



6(1 + go) + Jo -6 



kc^ 

kc^ 
-"o"o 



HnD 



+ 



(44) 



Combining these formulae we find that the luminosity distance as a function of D, the 
physical distance travelled is: 



dL{D) = D{i + 



3 [H.D 



+ 



kc^ 

ll + 4?o-T7272 



^o«5. 



mo 



1 

+ 24 



50 + 40go + 5jo - 10 



kc^ 
-"o«o. 



+ 



HnD 



(45) 



Now using the series expansion for for D{z) we finally derive the luminosity-distance 
version of the Hubble law: 



c z 



dL{z) = —{l+-[l 



Ho 



1 



<Io\z-- 



Qo - 3go + Jo + 



+ 



1 



2 - 2go - 15g^ - 15g^ + 5jo + lOgojo + sq + 



2k^{l + 3go) 



24 ^ 

+0{z')\. 



^0% 



(46) 



The first two terms above are Weinberg's version of the Hubble law. His equation 
(14.6.8). The third term is equivalent to that obtained by Chiba and Nakamura [8]. 
The fourth order term appears to be new, and (as expected) depends linearly on the 
snap. Prom the derivation above it is now clear that the 0{z^) term in this luminosity 
distance version of the Hubble law will also depend on the n-th time derivative of the 
scale factor. 

It is important to realise that this Hubble law, and indeed the entire discussion 
of this section, is completely model-independent — it assumes only that the geometry 
of the universe is well approximated by a FRW cosmology but does not invoke the 
Einstein field equations [Friedmann equation] or any particular matter model. Note 
that in comparison to the D{z) Hubble law, this di^z) Hubble law first differs in the 
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coefRcient of the 0{z'^) term — you will still get the same Hubble parameter, but if 
you are not sure which definition of "distance" you are using you may mis-estimate the 
higher-order coefficients (deceleration, jerk, and snap). The jerk jo hi'st shows up in the 
Hubble law at third order (order z^)\ but this was one of the parameters we needed to 
make the lowest-order estimate for the slope of the EOS. 

Warning: Not all authors use the same definition of the luminosity distance. In 
particular D'Inverno uses a definition that differs from Weinberg's by an extra factor of 
(1 -I- zY [18]. Weinberg's definition as presented above appears to the most standard, 
but if necessary the conversion is straightforward. □ 



5. A specific "a priori" model: Incoherent mixture of ly-matter 

Though the philosophy so far has been to avoid committing ourselves to any particular 
matter model, it is useful for comparison purposes to see how these ideas impact on the 
most popular models. A particularly common a priori model for the cosmological fiuid 
is an incoherent mixture of various forms of w-matter with each component satisfying 
a zero-offset equation of state: 

Pi ^vJiPi- (47) 

Integrating the conservation equation independently for each component of the mixture 
yields 

= po. (a/ao)-3(^+"'^) = po^ w, {a / a,)-^^'"-""^^ = l^o. {a / ao)-'^'+'"'\ (48) 

This model is sufficiently general to contain dust, radiation, cosmological constant, and 
standard forms of quintessence. Then 

Wo = ^ o — = w, (49) 

is simply the weighted average value of w. (A sum of w;.; over all matter components i, 
weighted by their present-day contribution to the Q parameter.) Similarly 

^ ^ dp/da|o ^ Y^i^oi Wi{l+Wi) ^ + w ^ aj ^ _ 
° dp/da|Q ^■iloi{l + Wi) l + uJ 1 + w 

That is, Wq (and hence go) provides information about the weighted average value of the 
Wi, while (and hence jo) provides information about how much spread there is in the 
various Wi. Uncertainties in the jerk parameter jo (which in a generic model manifest 
themselves as uncertainties in the slope parameter kq) in this specific model show up as 
difficulty in determining the weighting parameters Qqi- (For similar comments, see [12].) 
Note that the positivity of cr^ implies a constraint 

(l + w) (ko-w) > 0. (51) 
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If this inequality observationally fails, it means that the cosmological fluid cannot 
be described by any possible linear combination of tw-matter. In the approximation 
Hoao/c^ 1 this reduces to 

io>go(l-?o)- (52) 

At the next highest level in the Taylor expansion, again making use of the implicit 
function theorem, we have 

d^p p" - K p" 



dp2 (p')2 ■ 
After a brief computation this leads to 



d?p 
di^ 



w^{l -\-w)-\- w'^{l — w — w"^) — w 



2 



In terms of standard deviation and skewness. 



3 



(53) 



(54) 



Syj = {w — w)^ = — w + 2w^ — — 3(7^ — w^, (55) 

this yields 

d'^p 
dp 



PO J, 2 



S^il + w)-ai- 2al - 2al - Aa^ w - 2w' - 2w^ ^^^^ 



(l + tZJ)3 

Thus d^p/dp^lg is related to skewness in the distribution of the Wi. The general message 
to be extracted here is that the n-th Taylor coefficient in the EOS depends linearly on 
the (n+l)-th f2-weighted moment of the Wi. 

6. The observational situation 

As of March 2004 the observational situation is best summarized by the dataset 
described in the recent article by Riess et al. [19]. (See also the brief analysis in 
Caldwell and Kamionkowski [20]). The Riess et al. analysis presupposes a A: = 
spatially fiat universe (effectively, they adopt the inflationary paradigm and use the 
approximation c/{Hoao) <^ 1). They report that the jerk jo is positive at the 92% 
confidence level based on their "gold" dataset, and is positive at the 95% confidence 
level based on their "gold+silver" dataset. No exphcit upper bounds are given for 
the jerk, nor are any constraints placed on the snap Sq. The allowed region (for their 
preferred parameterization of the data in terms of [dg/d;2]o) is presented in their Figure 
5, and can usefully be summarized as follows: the allowed region is a narrow ellipse 
approximately centred along the line 

d-z = -'--2'' ^''^ 
and bounded by the boxes 

dq 



go e (-1.3, -0.2); 



G (-0.2, +3.8) (99% confidence, "gold"), (58) 
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and 



go e (-1.4, -0.3); 



G (+0.0, +3.9) (99% confidence, "gold+silver"). (59) 



To extract more detailed information, one needs to translate from [dg/d;2]o to jerk Jq. 
Working to third order in redshift it is straightforward to calculate {k — Q): 

A ( \ ( 1-qo 1 - + [dg/d^]o 2 . 3x1 /«n^ 
dL{z) = Jf^y + ^ g ^ +0{z)'>, (60) 



and so obtain 



2 ^9 
jo ^qo + 2qQ+ — 



(61) 





Using wq = —{1 — 2go)/3 it is then easy to determine 

Wo e (-1.2,-0.5); Jo e (-0.3, +5.9) (99% confidence, "gold"), (62) 
and 

Wo e (-1.3,-0.5); Jo e (-0.1, +6.4) (99% confidence, "gold+silver"). (63) 

Note that even at the 68% level both data sets exhibit a thin sliver of parameter space 
compatible with go < —1- (That is, Wo < — 1, so that the grand total of all stress-energy 
contributing to the cosmological fluid is "phantom" .) Also note that the upper bounds 
on the jerk are comparatively weak. 

Now consider the linearized equation of state using 

This parameter is very poorly constrained 

Ko e (-00, -5.5) U (-0.5, +00) (99% confidence, "gold"), (65) 

and 

K,o G (-00, -4.5) U (-0.5, +00) (99% confidence, "gold+silver"). (66) 

Note that very little of the real line is excluded, and that the entire positive axis is 
allowed by the data. Since this might at first be a bit surprising, let me explain the basic 
reason for this behaviour: If the permissible region contains a subset compatible with 
Jo 7^ 1 and go = —1, then this implies a subset of the permissible region is compatible 
with kq — ±00. That is, whenever there is a part of parameter space compatible with 
phantom matter, then kq ^ ±00 at the edge of the phantom region. 

Similarly, as long as the permissible region contains a subset compatible with jo = 1 
and go ^ —1, this implies a subset of the permissible region is compatible with kq — 0. 
Now jo = 1 is equivalent to 
dg 



dz 



= (l + go)(l-2go) (67) 



and it is easy to see that a subset of the Riess et al. permissible region is compatible 
with this constraint. With both of these conditions being satisfied the entire positive 
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real line is allowed for kq. The portion of the negative real line attached to — oo is 
associated with phantom matter, while the portion of the negative real line attached to 
the origin is associated with non- phantom matter (wq > — 1) with the peculiar property 
dp/dp < 0. 

It is also easy to see why at least some values of kq are excluded. If the permissible 
region contains the point — io = +1 (equivalently go = [dq'/d2;]o = 0) 
then, since this corresponds to Kq taking on the indeterminate value 0/0, it is easy to 
convince oneself that for any open set surrounding this point the computed values of 
Kq completely cover the entire real line — both positive and negative. But the point 
go = —1, Jo = +1 (^0 = —1, [dg/d2;]o = 0) is in fact excluded from both gold and 
gold+silver datasets at more than 99% confidence, which is ultimately the reason that 
at least some values of kq can be excluded. However, it should be emphasised that the 
constraints on kq are best described as extremely weak. 

A slightly different analysis (using the same raw data analyzed in somewhat different 
fashion) is presented by Caldwell and Kamionkowski [20]. Their result can usefully be 
summarized (again assuming A; = 0) as: 

go e (-1.1,-0.2); jo e (-0.5, +3.9) (95% confidence), (68) 

whence 

Wq e (-1.1, -0.5) (95% confidence), (69) 

and 

Ko G (-CX), -10) U (-0.6, +CX)) (95% confidence). (70) 

While the numbers are slightly different, the overall message is the same: the jerk is 
relatively weakly constrained and hq is very poorly constrained. 

If one wishes to constrain kq to at least lie in some bounded region of the real line, 
then one would need to seek improved data that might exclude the possibility qQ — —l 
{wo — —1). [This would then be the death-knell for phantom matter.] Similarly, if 
one wishes to constrain kq to at least lie in the "physically most reasonable" region 
(some bounded region of the positive real line), one would need to seek improved data 
that might exclude the possibility jo < 1- Obtaining a dataset of such quality would 
be extremely challenging: Assuming no change in the location of the center of the 
currently determined permissible region, this would correspond to contracting the 99% 
confidence intervals inwards to lie somewhere inside the current location of the 68% 
confidence intervals. This corresponds to approximately a three-fold decrease in the 
standard deviations obtained when fitting the Hubble law to the dataset; despite the 
advances in "precision cosmology" , this appears well beyond current (or even reasonably 
foreseeable) capabilities. 

Furthermore, one should note that the above considerations are subject to having 
additional data (or theoretical prejudices) available to justify setting k — [or more 
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precisely, c/{HQaQ) ^1]. In the absence of such data, the third order term in the 
luminosity distance Hubble law only implies [20] 

k 

go e (-1.1, -0.2); jo + e (-0.5, +3.9) (95% confidence). (71) 

Constraints on the equation of state, kq, are now even weaker. 

Note that appealing to the quartic term in the luminosity distance version of the 
Hubble law will not help, as that fourth order term brings in an additional free parameter 
(the snap Sq), so that there are still more free parameters than there are coefficients that 
can be measured. To side-step this particular problem the two possibilities are: (1) The 
traditional one, find additional data, above and beyond the Hubble law, that somehow 
constrain the space curvature k/al- (2) More challengingly, use a modified Hubble 
law that is independent of space curvature k/a^. For instance, the "physical distance 
travelled" version of the Hubble law D(z) presented in equation (34) is completely 
independent of space curvature — the challenge in this approach is to find a useful way 
of observationally measuring D. 

7. Discussion 

There are currently many different models for the cosmological fiuid under active 
consideration. Though these models often make dramatically differing predictions in 
the distant past {e.g., a "bounce") or future {e.g., a "big rip") there is considerable 
degeneracy among the models in that many physically quite different models are 
compatible with present day observations. To understand the origin of this degeneracy 
I have chosen to rephrase the question in terms of a phenomenological approach where 
cosmological observations are used to construct an "observed" equation of state. The 
key result is that even at the linearized level, determining the slope of the EOS requires 
information coming from the third order term in the Hubble law. Unfortunately, while 
the experimental determinations of the parameters appearing in the Hubble law are 
certainly improving, we are still somewhat limited in what we can say concerning the 
third-order term. Despite the fact that some parameters in cosmology are now known to 
high accuracy, other parameters can still only be crudely bounded [21]. In particular, the 
jerk is relatively poorly bounded, and as a consequence direct observational constraints 
on the cosmological EOS (in the form of measurements of Kq = [dp/dp]o) are currently 
extremely poor and are likely to remain poor for the foreseeable future. 
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